APPELL POLYNOMIALS AND THEIR ZERO ATTRACTORS 



ROBERT P. BOYER AND WILLIAM M. Y. GOH 



Abstract. A polynomial family {p n (x)} is Appell if it is given by = Y^m=o Pn{ x )t n or > ecm iv- 
alently, p' n (x) — p n -i{x). If g(t) is an entire function, g(0) 7^ 0, with at least one zero, the 
asymptotics of linearly scaled polynomials {p n (nx)} axe described by means of finitely zeros of 
g, including those of minimal modulus. As a consequence, we determine the limiting behavior of 
their zeros as well as their density. The techniques and results extend our earlier work on Euler 
polynomials. 
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1. Introduction 

Let g(t) be an entire function such that (7(0) 7^ 0. 

Definition 1. The Appell polynomials {p n (x)} associated with generating function g(t) are given 
by 

xt 00 

(i) ;br|> wt " 

Some important examples are: the Taylor polynomials of e x , with g(t) = 1 — t; the Euler 
polynomials, with g(t) = (e* + l)/2; and Bernoulli polynomials, with g(t) = (er — l)/t; and their 
higher order analogues. 

The asymptotics and limiting behavior of the zeros of these families have been investigated by 
many people; for example, [2], [7], and so on. 
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In this paper, we obtained the asymptotics and the limiting behavior of the zeros for all Appell 
families provided the generating function g(t) satisfies one further condition: that g must have at 
least one zero. We use the ideas in our earlier paper |2]; furthermore, we simultaneously simplify 
and generalize some of the techniques there. 

We found that the asymptotics in the general case are build from the basic example g{t) = 1 — t 
which coincides with the classical work of Szego on the Taylor polynomials of the exponential 
function. In our paper [2] , we found that the asymptotics for the Euler and Bernoulli polynomials 
are controlled by certain roots of g(t), the ones of minimal modulus. In the general situation, 
as expected, the minimal modulus roots of g(t) are needed to describe the asymptotics but there 
may be finitely many other roots needed to determine the asymptotics. These additional roots are 
determined through a geometric condition described in the terms of rotated and scaled versions of 
the Szego curve: |:re 1_:r | = 1, \x\ < 1, x £ C (see Figure [3]). 

We frequently use the following notations. Let Z{g) denote the set of all zeros of g and let 
i"o < r i < r 2 < ■ ■ ■ denote the distinct moduli of these zeros in increasing order. 

Recall that if K\ and K2 are two non-empty compact subsets of C, then their Hausdorff distance 
is the larger of sup{<i(x, K\) : x 6 K2} and sup{d(x, K2) : x £ K{\. 

Definition 2. For a family {q n (x)} of polynomials whose degrees are increasing to infinity, their 
zero attractor is the limit of their set of zeros Z{q n ) in the Hausdorff metric on the space of all 
non-empty compact subsets of the complex plane C. 

In the appendix, we discuss how the zero attractor is found in terms of the limsup and liminf of 
the zero sets. 




-0.4 -0.2 0.2 0.4 



Figure 1. Zeros for degree 1000 polynomial, with generating function g{t) = Jo{t) 



There is related work on the asymptotics and zeros of the Taylor polynomials for linear combi- 
nations of exponentials ^Cj-e AjX where the parameters Xj satisfy a geometric constraint pp. The 
techniques of proof are very different from our approach. 
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2. The Generalized Szego Approximations 

It is convenient to collect together several results from [2] and some extensions of them concerning 
the asymptotics of S n (x) = ^2k=o x The domains of where their asymptotics hold are critical 
in understanding the behavior for the Appell polynomials. 

Proposition 1. (Left-Half Plane) Let 1/3 < a < 1/2 and 1 < j. On any compact subset K 
of {w : $tw < 1}, we have 

V ; e nw v^(l-w) ' 

(2) Di-\w- n S n - X (nw)) = Di-\w- n e nz ) - ^=^7, g (l + (V" 3 ")), 

V27rn (1 - w) J 
where the big O constant holds uniformly for x £ K . 

The proof of part (a) is in [2] . Part (b) follows from an application of the saddle point method. 
The following Proposition is also from (2]: 

Proposition 2. (Outside Disk) Let S be a subset contained in \ w\ > 1 with distance 5 > from 
the unit circle, and let a be chosen so 1/3 < a < 1/2. Then 

e nw ^n~(w - 1) V V 11 ' 

where the big O term holds uniformly for w £ S. 

Proposition 3. (Evaluations of Integrals) If e < \w\ and j > 1, then we have 
1 f fe xt \ n 1 



2vri /| t | =e V f J t-w 



where D w is the differentiation operator 

Proof, (a) By expanding l/(t — z) into an infinite geometric series and performing a term-by-term 
integration, we obtain 



1 /■ /e xt \ n 1 



2m ./| t | =£ \ t / t - 2 



(it 




z2iri 

By the Cauchy integral theorem the terms correspond tom>n vanish. Hence 
1 /■ /e x *\ n 1 



dt. 



- E - - 

7 £—1 7 m \ 2-7T 



dt = > r 4 e zm t -n+m rft 



27ri \ t I t — z z ^— ' z m \ 27U /ui_- 

l'|- £ v 7 n-l>m>0 V -'Itl-e 



1 1 ( xn ) 



E 



n—m—l 



n-l>m>0 y ' 

-\ z -n+i V" {xnz) 



1 E 



n—m—l 



z ^ (n — m — l)! 

n-l>m>0 v ' 

_ z _ n ^ (xnzp = { _ l)z -n Sn _ i{zxn y 

n-l>j>0 



Part (b) follows from differentiating (a) j — 1 times with respect to z. □ 

3. Asymptotics Outside the Disk L>(0;l/r ) 
Theorem 1. Let K be any compact subset in the annulus A{\/rQ]oo). We have 

1 (l + (l/„)). 
(xe) n /V2im g(l/x) V y 1 " 

holds uniformly for x £ K. 

Proof. We shall find an asymptotic approximation to p n {nx) in the region A(l/ro; oo) = |x : \x\ > i|. 
Use the generating relation equation ([I]) to get 

1 I e 3 * 

Vn\X) = — ^ rr «ii 

where e < ro- Since both sides of the above equation are entire functions of x, by analytic con- 
tinuation this representation for p n {x) is valid for all x 6 C. Hence we can replace x by nx to 
get 

The above expression is valid for < e < ro and is the starting point of the analysis in the sequel. 

Let K be an arbitrary compact subset C {x : \x\ > — } and let x £ K. We can certainly choose 
e small enough so that for all x G K, \ex\ < 1. By a change of variables, we get 

. , x n f /e 4 \ n dt 
p n (nx) = <b — — ■ 

2^ J\t\=e\x\ V * / 

Observe that the zeros of g{t/x) have the form ax where a G Z{g). Moreover, they must lie outside 
the closed unit disk since \x\ > 1/ro, so we can deform the integration path from the circle with 
radius e\x\ to the unit circumference. Thus 

x n f fe l \ n dt 

p n {nx) ■- 



(2) p n (nx) 



2tti yj t | =1 V t J tg(t/x 
%n f n(t-Int) dt 



2ni 7J t | =1 tg(t/x) 

It can be easily seen that t = 1 is the saddle point of the integral and the classical saddle point 
method is applicable here [3]. Hence 

where the implied O constant holds uniformly for x £ K. □ 
The last equation can be written as 

1 + O(-) , \x\> l/r . 



We have the: 

Corollary 1. (a) On the complement of the disk D(0; 1/ro), nm — hi 

n— >oo 77, 

i/ie limii ZioWs uniformly on compact subsets. 

(b) The zero attractor must be contained in the closed disk D(0; 1/ro). 
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Pnjnx) 

nn 



where 



Note that part (b) follows easily from (a) since g{x) never vanishes on the disk D(0; 1/ro). 



4. ASYMPTOTICS ON THE BASIC REGIONS Re 

Let ro,ri,. . . denote the distinct moduli of the zeros a of the generating function g. Fix an 
integer I. We fix p > so it is not equal to any zero modulus {ro, n, . . . }. For each zero a G Z(g) 
with \a\ = re, we consider the circle \x\ = l/\a\ and the disk D(l/a,S a ). 

Now the tangent line T a to the circle |x| = l/\a\ at the point 1/a determines the half-plane H a , 
which contains 0; that is, 3?(ax) < 1. We choose ee > to be less than the distance from the 
portion of the tangent line T a that lies outside the disk D(l/a; d a ) to the circle |x| = l/\a\ for any 

\a\ = r i+1 ; that is, ee < \Jl/rf + 5% - l/r £ . Finally, we make the requirement the disks D(\Ja\ 5 a ) 

be mutually disjoint for all a G Z{g) with \a\ < p. 

Definition 3. With these conventions, the region Re is described in terms of the half-planes H a 
and disks as 

(3) Ri = f) {H a \ D{\- S a ) : \a\ = n) \ D(0; ^ + 

We note that the regions Ri are not disjoint; in fact, by construction, its inner boundary which 
consists of the portion of the circle \x\ = ^— j- + e that lie outside the disks D(l/a; 5 a ), \a\ = ri + \, 
actually lies inside the region Re + \. 

It is convenient to introduce a region that contains all of the Rps. 

Definition 4. Let R p be the domain given as 

(4) R p = p| {H a : a G Z(g), \a\ = r } \ [(J {D(l/a; 5 a ) : a £ Z(g), \a\ < p] U D(0; l/p) . 

Note the order of dependence: first we can given the cut-off modulus p > for the moduli of the 
zeros; next, S a > for each o G Z{g) is given and is a function of p [see later section], then finally, 
ee is determined relative to each zero moduli ri which is a function of 5 a . 

For any a G Z(g) with ro < \a\ < p, let s a (t) be the singular part of 

1 

tgjt) 

at its pole a. Next we set gi(t) to be 

(5) gi(t) = - : a G Z (9),r < \a\ < p} 

we see that g±(t) is analytic for |i| < p. 

We develop the asymptotics for {p n (nx)} on the regions Re where ro < re < p. Now we saw 
already that we can write p n {nx) as 

where s(t) = ^{s a (t) : a G Z(g),r < \a\ < p}. 

Lemma 1. With gi(t) given above in equation^ we have 

i [■ / p xt\ n r n-l„n 

M / ( — ) 9i(t)dt= ^(l/ar) (1 + 0(l/n)) 

uniformly on compact subsets of the annulus ^4(l/p;oo). 
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Proof. Let x £ K C A(l/p, oo). By a change of variables, we write 
1 f fe xt \ n x n ~ l f I V\ n 

By construction, g(t/x) is analytic on a disk of radius greater than 1. So the contour in the last 
integral can be deformed to the unit circle |x| = 1 without changing its value. Finally, by an 
application of the saddle point method we find that 

„n— 1 r / p t\ n ™n— l„n 

*«L(V «.(«/.)* -^«.d/.)(i + o(i)). 

□ 




Figure 2 . Generic Plot of Polynomial Zeros and Zero Attractor When g Has Two 
Roots; Tangent Lines and Circles Displayed 



To state the next two lemmas, we need to introduce special polynomials I n (z) in z^ 1 and J (a; z) 
in z. 

The polynomial I n (z) comes from expanding the derivative of D" 1 ~ 1 (z~ n e nz ). Consider 

m ~ 1 / — 1\ 

D^- 1 {z- n e nz ) = Y.( m , ){Dlz- n )(D n z l - l - p e nz ) 



m—l 

E 

p=0 



m — l 



P 



{-n){-n - 1) • • • (-n - p + l) Z - n - p (n m - l - p e nz ) 



m—l 



z -n e nz n m- 



~ l ( m 1 ) (-")(-" ~ 1) ■ • • (~n- P + l)(nz)- p 



nz) 



z -n e nz n m-lj m _ i{nz)i 



where I m -i{z) is given in 
Definition 5. 

(6 ) w*)-£(-w( ra ; 1 )( n+ *- 1 )^. 

For a G Z(g), we define J(a; z) which are also polynomials in We write out the singular part 
s a (t) of the function at its nonzero pole a by 

(7) s a {t) := V 

m=l 7 

where /3 a is the order of a as a zero of g(t) so b a ^ a ^ 0. 

Definition 6. For a S Z{g), let J(a;z) be the polynomial in z given as 

(8) J(a;z) = £-^^ 1 / m _ 1 (az). 

m=l 

Lemma 2. Let a 6 Z(g) and let x £ K, a compact subset of the half-plane H a , $t(ax) < 1. Then 

^- / V s «(*) * = -«~ n e nax J(a; nx) + ^= Sa (l/x) (l + OCn 1 - 3 ")) 
27H 7| t|=e V f / V2vrn 

where s a (t) is the singular part of l/(tg(t)) at the zero a of g(t). 

Proof. We first write out the integral in terms of the singular part s a (t) 

I I' 

1m 



where the coefficients b ajTn are given in equation ([7j). We now study the asymptotics of the typical 
term D^- 1 (a- n S n - 1 (nax)). 

We may use the generalized half-plane Szego asymptotics with | < a < \ because of the 
restriction that a £ Z{g) with \a\ < m to obtain 

D™- l {or n S n ^{nax)) = x n+m ~ x D^~ x {{ax)' n S n - X {nax)) 



x n+m—l 



Combining these estimates we obtain 

D^ia^Sn^inax)) = x n+m - 1 {(ax)- n e nax n m - 1 I m ^ 1 (nax) 



y/2m~i (1 - ax) 1 



1 + 0{n l ~ Za ))} 



(m — 1 V P n T n+m—l 

(9) = a-»e""(nx) m - 1 / m _ 1 (na aJ ) - (l + O^ 3 ")) 

y27m (i- — ax) 



Hence after summation we obtain 

^xt \ n A» 



1 f / Xt \ n " a L 

f (t) Sait)dt = -Era^VVjH) 



2vri 



n n-l i 

(10) = - a -V^J(a;nx) + ^=s (-) (l + Ofa 1 - 30 )). 

V27rn x 



□ 



Corollary 2. For a G Z(g), \a\ < re, we have 

1 P / xt \ ^ Tl 71 — 1 

ir-. I [—) s a (t) dt = -a- n e nax J{a- nx) + ^= S(1 (l/x) (l + 0(n^ 3a )) 
27rz 7| t | =e V t J \/2irn 

uniformly on the compact subsets of Re, where s a (t) is the singular part of l/(tg(t)) at the zero a 
ofg(t). ' 

Lemma 3. Let a G Z(g) and let x G K, where K compact subset of the disk- complement 
A(l/|a|;oo). Then 

1 f / ' p xt \ n p n T n-l 

^L(t) «.W* = ^Mi/.)(i + o(„— )). 

Proof. We will use the disk-complement generalized Szego asymptotics. For z in the annulus 
A(l + c, oo), for any c > 0, we have 

z n r e n(C-lnC) 

S n --v{nz) = -— & — d( 

2tti J\q =1 (-Z 

By Dividing z n and taking derivatives up to order m — 1, we get 

f m _ i V r P "(C-inC) 

2 V U " 2vri jT[ C | =1 (C-^) m 

In the above, replace z by ax to obtain 

D™-\ a -"S n - 1 (nax)) = *"+™- 1 ZC- 1 ((a*)-"S„_ 1 (™>*)) 

(11) = (l + (»'^)). 

V 7 (1 -ax) m v v 

By summation, we obtain the asymptotics for the original integral: 

. — I sJt) dt = - Y 

e n x n - 1 g 1 
\f2wn a x 



1 /* / Xt \ n $ a 7, 

n n-l i 

(12) = ^^^(^(l + O^- 3 ")). 



□ 



Corollary 3. For a G iwit/i r^ + i < |a| < p, we /tave 
1 f f p xt \ n p n n-l 

V S ^ dt = ^|=*a(l/x) (1 + 0(n^)) , 
2vr? 7| t | =e V * / V2vrn 

uniformly on the compact subsets of Re- 

The remaining case for the above integration involving s a (t) on the disk D(l/a; 5) will be handled 
in a later section. 

Proposition 4. For x e Re, we have 
p n (nx) 



= -, — -— — V2nn { J(a; nx)— — — : a G Z(q), \a\ < 

(xe) n /V2^i x gi (l/x) ^\ V V(oi)« ~~ 

+ E { ^ Q) : a £ |a| < p} (1 + 0{n l ~^)) 



re 



uniformly on the compact subsets of Rg, where <fi(x) = xe 1 x and 1/3 < a < 1/2. 

Proof. Putting the last two corollaries into Equation ^ and simplifying, we have 

/ n x n ~ l e n .1. / ^.1. 
Pn (nx) = —== 9l {-) 1 + O - 
V27rn x \ n 

~ Yl {a~ n e nax J(a; nx) : a £ Z(g), \a\ < p) (l + 0{n l ~^ a )) . 



Proposition 5. For x G Rg, we have 

Pn{nx) = ]L _ nx )^_ . a g Z ( ) | a | < } + 0( n l-3a) 

uniformly on the compact subsets of Ri, where <f>(x) = xe 1 ~ x and 1/3 < a < 1/2. 
Proof. By the definition of g\{t) (see equation ^5|), we see that 
1/1 

13 -gi - 

x \x 



)=g^)-^{ 1 x Sa{1 x ):aeZ{9)M<P }- 



Insert this into the above. Since the s a {\/x) term cancels, we have uniformly for x € Rf 
p n (nx) 1 



(14) 



(ex) n /\f2im g(l/x) 

-V2^Y1 {{axe l ~ ax )- n J{a; nx) : a G Z(g), \a\ < p} + O (ra 1 " 3 ") . 



Lemma 4. If a G Z'(fl') mt/i |a| < p and x £ Rg, then 

J( ° ;ra) = (S^ (ra) ^ 1 (^ i ) & " (1 + ° (1)) ^ 

Proof. Recall that 

J(a;nx) _ * W-t^) , Jm _ l(nol) = £)(-!)-('» - *) + 

m=l v P=0 \ ^ / \ ^ / 

It is easy to see that 

'n + p — 1\ , (ax) 



p J ' p! 



(l + o(l)), 



that is, as n —* oo 

m— 1 / 1 \ / 1 \ m— 1 



□ 



□ 



p=0 

Hence 

(15) J{a,nx) = T p^{nx)^( a ^) fia \l + o(l)). 

[Pa ~ 1)! \ ax J 

Since the coefficient b a ^ a in the definition of the singular part s a (t) is nonzero, we find for fixed x 
that the precise order of J(a;nx) as a polynomial in n is n^ a_1 . □ 

We note the following 



Corollary 4. lim — = — — — , x G Re provided \<p(ax)\ > 1 for all a G Z(g) imi/t 

(xe) n /V27rn 5(l/a0 

|a| < 77. 

p n (nx) 



— ln|0(aox)|, /or a// x G Re satisfying \cf)(aox)\ > 1 for 



Corollary 5. lim —In 

rwoo n (xe) n /V27rn 

ao G ^(fl 1 ) 1 4>o{ ax )\ < l</>( aa; )l / or a M a a o such that \a\ < re- 

By construction, Re C R p for all £ chosen so re < p. Consequently, we have the asymptotics: 

Theorem 2. On R p , we have the following uniform asymptotics 

Pn{nx) 1 

; = — ; — rr 1 + ^ 

(ex)«/V2^ 9(l/x) 

-y/2Kn^{(axe l ~ ax y n J(a-nx) : a G Z( ff ), |a| < p} + O (n 1 " 3 ") , 
w/iere 1/2 < a < 1/3. 

It remains to develop the asymptotics in the disks D(l/a; 5 a ) and well as determining domination 
among a G Z(g) of |(^>(ax)|. 

5. Dominant Zeros and Szego Curves 

Let (f>(x) = xe 1 ~ x which is an entire function that is conformal on the unit disk. The standard 
Szego curve S is the portion of the level curve \<j>(x)\ = 1 that lies inside the closed unit disk or 
equivalently, inside the closed left-hand plane 5ft(x) < 1. S is a closed simple closed convex curve; 
in fact, it has the form t = ±\/ e 2 ( x ~ s ) — s 2 where x = s + it and s G [— W(e~ 1 ), 1] and W is the 
principal branch of the Lambert VF-function. 

Definition 7. Let a be a nonzero complex number. We call any curve of the form a Szego 
curve. 




Figure 3. Szego Curve: \ze 1 ' z \ = 1 and \z\ < 1 
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Lemma 5. Let a,b be non-zero distinct complex numbers. Then the intersection -S D tS has at 
most 2 points. 

Proof. The intersection of the two curves \S n \S must satisfy |^(ax)| = |0(6x)|. This modulus 
condition determines a line; so, the intersection of the two curves lie on the and a line. Since 
^5 is convex, the intersection contains at most two points. 

Write x = s + it and b — a = a + i(3. Then |<^>(aa;)| = |0(6a;)| reduces to the line: 

\ae- ax \ = \be- bx \, \e^ x \ = \b/a\, \e^ ib - a)x] \ = \b/a\, 
— a)x] = In \b/a\, as — (3t = In \b/a\. 

□ 

Definition 8. Call a E Z{g) a dominant zero if either \a\ = ro orl/a does not lie in the interior 
of any of Szego curves where b is a minimal modulus zero, \b\ = tq. IfX/a lies on one of these 
Szego curves, call a an improper dominant zero; otherwise, a is proper. 

Note: if a is an improper dominant zero, then there exists a minimal modulus zero b so the 
intersection n ^5 consists of a single point, namely l/o. In [1, the behavior of the zeros of the 
Taylor polynomials of Sj=i Cje^i z is determined if the numbers Xj satisfy a condition in terms of 
a convex m-gon. For the Appell polynomials, the geometric condition is more subtle. 

Let W denote the principal value of the Lambert W-function. Since the radius of the largest 
circle centered that lies in the interior of the standard Szego curve S is W(e~ l ), we have the 
following 

Lemma 6. If a' is a zero of g such that \a'\ > ro/W(e^ 1 ), then a' must be a non-dominant zero. 
Hence, there are at most finitely many dominant zeros. 

Definition 9. // a and b are two dominant zeros such that n consists of two points, then 
the line \4>(ax)\ = \(p(bx)\ determines two half-planes. Let E a + b denote the closed half-plane that 
contains 1/a. 

Let ai,a2, ■ ■ ■ ,a n be the dominant zeros of g. Recall that if f(x) is any analytic function on a 
domain D then its modulus is a subharmonic function on D. We need two basic properties 

of subharmonic functions: they satisfy the maximum modulus principle; and the maximum of two 
subharmonic functions is still subharmonic. 

Definition 10. Set <&(x) = max(|(/>(aix)| _1 , • • • , |(/>(a n x)| _1 ) so is subharmonic on the punc- 
tured complex plane. 

We work with |</>(aa;)| _1 rather than |<^(aa;)| so we may apply the Maximum Principle since 
|^(ax)| vanishes at x = 0. 

We observe that the level curve = 1 divides the complex plane into finitely many connected 
components just as the original curve |0(x)| = 1 divides the complex plane into three connected 
components. We make the: 

Definition 11. Consider the connected components of the complement of the level curve = 1. 
Let Dq be the closure of the connected component that contains 0. Note that if x € Dq and is x ^ ; 
then <£(x) > 1 with strict inequality when x lies in the interior of Dq. 

We assume that p > is chosen so large that the closed disk D(0; 1/p) is a subset of the interior 
of Dq. Furthermore, the singularity in log |0(ax)| — log |</>(6x)| always cancels so this difference is 
always a harmonic function on C. 

Let a' be a non-dominant zero of g such that \a'\ < p. Then we know that -^<S is a subset of the 
interior of Do and that the disk D(0; 1/p) lies inside of ^<S. 

The following Proposition follows easily from the definition of Dq and dominant zeros: 

li 



Proposition 6. Uniformly on the compact subsets of R p \ Dq, we have 

p n {nx) 1 

urn = — — — — 

(ex) n /v2vrn g{l/X) 

To understand the asymptotics inside -Do requires subharmonic function theory. 

By construction, this means that \4>(a'x)~ l \ < 1 for x G 9Dq, while $(x) = 1 for x G <9-Do- 
Before we can apply the Maximum Modulus Principle for subharmonic functions, we need to deal 
with the common singularity at 0. However, this is easily dealt with by multiplying both and 
by \x\ which shows that we can remove this singularity relative to the inequality. 

In particular, <£(x) is strictly larger than |</> _1 (x)| for all x G Dq, with i / 0. So there exists a 
positive constant, say a, so that 

> ^(a'x)!" 1 + a, x£D \ D(0;l/p). 

We state the above discussion formally: 

Proposition 7. Let a' be a non-dominant zero of g. On the domain Dq \ D(0; 1/p), we have the 
order estimates 

\4>(a'x)\- n = 0((*(x) - a) n ) = o($(x) n ). 

Definition 12. Given a proper dominant zero a of g, let 

D a = {x £ Dq : |0(ax)| < |0(6x)|, for all dominant zeros b ^ a, \(f>(ax)\ < 1}. 

Note that the definition of D a is independent of the choice of p sufficiently large. Further, D a 
has the alternate description in terms of the half-planes E a + 

— Interior(5) fl f^|{-E a + 6 : b is a proper dominant zero of g, b ^ a}. 

It is easy to state formally the basic properties of D a : 

Lemma 7. Let a be a proper dominant zero of g. Then D a is a non-empty compact connected 
subset of Dq. 

Proof. Since ^5 is a convex curve and the intersection of half-planes is connected, the set D a must 
be a connected convex set. □ 

We now restate Proposition [7] relative to a domain D a : 

Proposition 8. Let a' be a non-dominant zero of g with \a\ < p. Then a' must lie in a domain 
D a for some proper dominant zero a of g. For x G D a \ D(0; 1/p), the following holds uniformly 

\4>(a'x)\- n = 0{\(j){ax) - a\~ n ) = o(\4>{ax)\- n ) . 

Finally, we now have a refinement of Theorem [2] as 

Theorem 3. Let p be chosen greater than l/\a\ where a is any proper dominant zero of g. Then 
on R p , we have the following uniform asymptotics 

(ex) n /V2^i g(l/x) V y ' " 

nx) : a G Z(g) and dominant } + O (n l 3a ) + o($(x)), 

where 1/3 < a < 1/2 and = max{|0(ax)| _1 : a G Z(g) and dominant }. 
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6. Asymptotics Inside the Disk D(l/a';5 a ') When a' is a Non-Dominant Zero 
We first state an easy consequence of a previous Proposition. 

Proposition 9. Let a' G Z(g). Then on the disk D(l/a'; 5 a /), the normalized polynomials have the 
asymptotics 

; = — 9i — 1 + ^ — + > w(a; nx )—< — r~ : a £ Z(q),\a\ < \a \ ,a ^ a > 

+ E {^ Sa Q) : a G Z (5 f )> a / |a| < p| +v a >(x), 

where 

Pa' 



m=l 



Proposition 10. Let a' 6e a non-dominant zero of g with \a'\ < p. Then there exists a choice of 
5 a r > such that 



where p is the cut-off bound for the zeros of g. 

Proof. To estimate a a , we make use of the elementary estimate: If f(z) is analytic function of z, 
then for any e > 0, we have 

\D{^f{z)\< { ^— ^ max |/(C)| . 

e J \C-z\=e 



By the definition of a a > (x) , we find 



K'l < 



Pa' . 

On/, 



< il^d ..max |C^n-i«*)| 

m=l <V IC-a'|=^ 

- ^a' , (\C\~ n S„-i(\Cx\n)) 

where ifj , > is a constant that depends on the zero a' and the radius 5 a i. 
To go further we observe for x G <5 a /) and |C — a'| = 6 a r. 

\(x\ < (la'l + <5 a ') |x| < |a'| \x\ + \x\ S a i 

< 1 + \a'\ 8 a i + |x| 5 a > = 1 + <5 a /(|a'| + 

Since |a'| < p by assumption, < 1 + 2p5 a '. But |£| > |a'| — S a i and |x| > ^ — <5 a /, so we get 

\(X\ > (|a'| - <y„0 - <J ^ > 1 - *a' + |«'|) > 1 - 2<W- 

Collecting these two inequalities, we get 

l-25 a 'p< \(x\ <l + 25 a 'p. 
13 



Now use that \S n -i(nt)\ < e nt : 

max \e(x\- n 5 n _i(|Cx|n) < e~ n \l - 25 a > p\~ n e< l+2 ^' p) 

K-a'\=6 a < 

= \l-25 a ,p\- n e 2nS «' p 

For < x < 1/2, 1/(1 — x) < e 2x ; if we choose 5 a > such that 25 a >p < 1/2, then we have 
|1 - 2S a >p\~ n < e 4S -' p . With this choice of 5, we obtain the desired bound 

max (|eCx|- n S n ^(\Cx\ n )) < e A5 »' p e 2n5 «' p = e 6n5 «' p 
\C-a'\=S a i 

□ 

7. Zero Attractor and the Density of the Zeros 

In our paper [2] , we determined the limit points of the zeros of the Euler polynomials by means 
of the asymptotics and the zero density. Here, we separate out first the question of find the support 
of the zero density measure, which is, of course, the zero attractor. Then we determine the zero 
density by applying our general result in the appendix. 

Proposition 11. Let f n (x) = \/2vrn p n (nx) j (xe) n . Then the following limits hold uniformly on 
compact subsets of the indicated domains: 

(1) On the domain A(l/ro;oo), lim — ln[/ n (x)] = 0. 

n— »oo TL 

(2) On the domain R p , lim — ln[/ n (x)] = 0. 

n— >oo TL 

(3) On the domain D a n A(l/p;oo) where a is any dominant zero of g, lim —ln[f n (x)] = 

n— >oo 77, 

— In (j){ax). 

Proof. We use the asymptotic expansions for p n (nx) developed in the previous sections. For |x| > 
1/ro, we noted already that the indicated limit must be 0. 

We observe that if a' is a nondominant zero of g with \a'\ < p, then for 6 a i > sufficiently small, 
the disk D(l/a!;8 a i) will lie in the domain D a for some dominant zero a, then on D(l/a';8 a i), 

lim — ln[/ n (x)l = — hx(f)(ax). □ 

n^oo n 

To describe the zero attractor requires a closer examination of the boundary of each domain D a 
where a is a proper dominant zero. 

The boundary dD a where a is a proper dominant zero of g has several natural families: dD a ndDo 
which is an "outer boundary" and a polygonal curve consisting of the line segments dD a n dD^ 
where b is another dominant zero of g. Note that dD a n dDf, is a subset of Dq. It will be useful 
to subdivide dD a n 8Dq into two connected components denoted by dD^ that come from deleting 
{1/a} from [dD a ndD Q ]. 

Lemma 8. The zero attractor of the Appell polynomials {p n ( n x)} must lie inside the compact set 

{dD a : a is a proper dominant zero of g} . 

Proof. First, we let x* let in th.6 infinite exterior of Dq- Recall that lim ri >OQ \/2irnp n (nx)/(xe) n = 

l/g(l/x) uniformly on compact subsets. If x nk is a zero of p nk (n^x) and x nk — > x*, then appealing 
to this limit we find that the limit must be while the right-hand side is l/g(l/x*) ^ 0. Secondly, 
suppose x* lies in the interior of Dq but not on any boundary set dD a , where a is a dominant zero. 
By construction, x will lie in the interior of one of the domains where b is a dominant zero. 
Then lirrin^oo | \j2itnp n {nx) j (xe)™| 1//n = |</>(6x)| uniformly on compacta in the interior of D a . By 
the same reasoning as before, x* cannot be a limit of zeros. □ 
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The following Theorem is an immediate consequence of the above lemma together with the result 
of Sokal in section A. 4 of the Appendix. 

Theorem 4. The zero attractor of the Appell polynomials {p n (nx)} is given by 

(^J {dD a : a is a proper dominant zero of g} . 

Proof. Let a be any proper ominant zero of g and let x* G dD^- Let e > be given. Then we find 
that 



lim In 

n—>oo 



p n (nx) 



(xe) n /\/2vrn 



lnl 



0, xeD(x*;e)\D , 
(ax) |, x e D(x*;e)Dlnt(D ) 



holds uniformly on compact subsets. Next suppose that x* is nonzero and lies on one of the line 
segments of the form dD a n dD^ where b is another proper dominant zero. Again, we find that 



lim In 

oo 



p n (nx) 



(xe 



) n /\ / 2vrn 



■]n\<j)(ax)\, x G D(x*; e) D Int(-D a ), 
-ln\<f>lbx)\, x e D(x*;e)f]lnt(D b ) 



which also holds uniformly on compact subsets. By Sokal 's result [B] which is described in the 
appendix, we conclude that x* is in limsupZ(p n ) since there can be no harmonic function v{x) on 
the disk D(x*; e) that satisfies the inequalities 



lim inf In 

n— >oo 



p n (nx) 



(xe) n /^/2vfn 



< v(x) < lim sup In 



p n {nx) 



(xe) n /\/2irn 



On the disks D(l/a';5 a i) where a' is a non-dominant zero which must lie inside Dq, the contri- 
bution of the dominant zero dominants. 

This reasoning handles all but finitely many points: 1/a where a is a dominant zero of g. However, 
since the zero attractor must be a compact set and points in D(l/a; e) n [dD a n 8Dq\ lie in the zero 
attractor, we conclude that 1/a also lie in the attractor. □ 

Theorem 5. Let a be a proper dominant zero of g. 

(a) The zero density measure on any proper subcurve of dD a ndDo is the pull-back of the normalized 
Lebesgue measure on the unit circle under the conformal map <p(ax). 

(b) Let b be a proper dominant zero of g so b / a. Then the zero density measure on any proper 
line segment of dD a n dDb is Lebesgue measure. 

Proof. For both parts, we can use the asymptotics given in Theorem [3j 

For part (a), let f n {x) = ^/2itng{\ / x)p n {nx) / (xe) n . Let a be a dominant zero of g, and let 
C be a proper subcurve of 8Dq n dD^. Then there exists a neighborhood U of C such that 
U C R p fl [(C \ -Do) U D a ] so that the asymptotics in Theorem [3] can be written as 

, ^TL- = -err-, (! + °(Vn)) - v& + 0{n^) + (*? >B (z)), 

(xe) n /v27rn giX/x) <p(ax) n 

where $i j0 ,(x) = max{l, |^>(ox)|}. Hence, by dividing by g(l/x)), we find that f n (x) has the form: 

f n (x) = 1 + a n (x)(f)(axy n + e n (x), a n (x) = -g(l/x)J(a; nx), 

where 

e n (x) = 



oil) 



x e UD{C\D ), 
o{(f>{ax)~ n ), x G U D (D a n R p ). 

Since (j)(ax) is conformal in the disk D(0;l/|a|), we may apply Theorem [6] from the Appendix 
Section [A.4| on the density of zeros. 
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Let a and b be two distinct proper dominant zeros of g such that dD a n dD b is nonempty. On 
D a Dbf] R p , the asymptotics in Theorem [3] can be written as 

Pn (nx) = + 0(l/n)) -y/2^(j( a;nx)—^— + J{b ;nx)—^— 

+ J(a'; nx) —r-r^r~ '■ a ' proper dominant zero, a' ^ a, b } ) + 0(n 1 ~ 3a ) + o(&(x) n ) 
<p(a x) n J 

= SIM (1 + 0(1/n)) - ^ ( J(a; + Jlk ~>i(£)0 + 0( "" 3 " ) + "WW*"' 

where ^^(x) = max{l/|0(ax)|, l/|^(6x)|}. 

Let L be a proper line segment of the intersection dD a n 9.D&. Let f7 be a neighborhood of L so 
both |0(ax)| < 1 and |0(foc)| < 1 for x £ U. On the intersection U n R p , we work with a different 
normalization than before: 

rp ( \ <p{ax) n 

J-n\X) = ~j= p n {nx). 

\/2nn(xe) n J(a; nx) 

Note that in this normalization the term that contains 4>(ax)~ n becomes the constant 1 for T n {x). 
Of course, this new normalization has exactly the same zeros as p n {nx) in U so the zero density is 
unchanged. Then we find that 

T n {x) = 1 + a n {x)i){x) n + e n (x), 

where 

4>{bx) b J (a; nx) 

and 

en(x) = - J^f - (OCn 1 - 3 -) + (*^(x))) • 
V 2TrnJ[a; nx) 

On J7, we have that ^(ax) n $" b (x) = max{l, \i()(x)\ n }; while on D a n U, \tp(x)\ < 1 and on D b n C/, 
|^>(a;)| > 1. This allows us to write e n (x) as 



e n (x) 



o(ip(x) n ), xeD a nu, 
o(i), xeD b nu. 



By construction, cj){ax) / (f)(bx) = |e^ b a ) a is a conformal map on f/ n i?p that maps L onto an arc 
of the unit circle. By Corollary [7] in the Appendix section A. 4 the result follows. □ 

We close with several examples that illustrate the main constructions in the paper. 

Example 1. Let g(t) be an entire function whose minimal modulus zero a\ = 1 such that all its 
other zeros a satisfy l/\a\ < W(e^ 1 ) ~ 0.27846. Then the zero attractor for the associated Appell 
polynomials coincide with the classical Szego curve in Figure |3j 

Example 2. The higher order Euler polynomials En\x), where m G Z + , have generating function 

g(t) = (e* + l) m /2 m ; while the higher order Bernoulli polynomials Bn(t) have generating function 
g(t) = (e — l) m /t m . Then their zero attractors are independent of m and coincide with a scaled 
version of the zero attractor for the Taylor polynomials for cos(x), see Figure |4j 

Example 3. The zero attractor for the Appell polynomials associated with generating function 
g{t) = Jo(t), where Jo(t) is the zero-th order Bessel function, is a scaled version as the zero 
attractor for the Taylor polynomials for cosh(x), see Figure [TJ since the minimal modulus zeros of 
Jo(t), a = 2.404825558, are the only dominant zeros since all the zeros of Jo(t) lie on the real axis. 
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Figure 4. Zero Attractor for Taylor polynomials of cos(x) 
Example 4. Let g(t) = (t — 1) {t 2 + 2). See FigurelHlfor its zero attractor and zeros for degree 400. 




Figure 5. (a) Zeros for degree 400 polynomial with generating function g(t) = 
(t — 1) (i 2 + 2); (b) Zero Attractor with polynomial zeros 



Example 5. Consider the Appell polynomials with generating function g(t) = (t — a)(t— b)(t — c) 
with a = 1.2e i3?r//16 , b = 1.3e l7?r//16 , and c = 1.5. In this case, all three roots of g(t) are dominant. 
See Figures [6] and [7j 

These last two examples both illustrate the following general fact. We assume that the generating 
function g(t) has exactly three proper dominant zeros a, b, and c. Then the three lines determined 
by |(/>(ax)| = |0(6x)|, |0(ax)| = |^(crc)|, and |0(6x)| = |^(cx)| have a common intersection point, a 
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FIGURE 6. (a) Zero Attractor only, for generating function g(t) = (t — a)(t — b)(t — c), 
a = 1.2e i3n / 16 , b = 1.3e i77r / 16 , c = 1.5; (b) Boundary of the Domain D . 
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Figure 7. Zeros for degree 400 polynomial together with the Zero Attractor, for 
generating function g(t) = (t - a)(t - b)(t - c), a = 1.2e i37r/16 , b = 1.3e i77r/16 , c = 1.5 

so-called "triple point." This follows by interpreting the lines as the boundary between the change 
of asymptotics of the Appell polynomial family; that is, the boundaries of the domains D a , Db, and 
D c . 

Appendix A. Density of Zeros 

A.l. Introduction. We generalize the density result for the zeros of the Euler polynomials in [2] 
to highlight how the asymptotic structure of the polynomial family may determine the density of 
its zeros. 
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Let ip(x) be an analytic function on a domain DcC that is conformal on D. We write £ = ip(x). 
We sometimes write x(Q for x = ^ _1 (C)- 

We assume that there exists eo > and < a < (3 < 2ir so that the annular sector 

(16) S = {pe w :p£ [l- eo ,l + eo ],0 G [ a , p]} 
lies in the image ip(D). Next we define two subsectors of S as 

S + = {pe* e : P £ [l-e o ,l),0€ [a, [3]} 

5_ = {peV-.pe (1,1 + eo],9€[a,0\}. 

Let C be the unimodular curve ip~ l ({e ie : 9 G [a, /?]}), so \(f>(x)\ = 1 for x G C. By construction, C is 
smoothly parametrized as x(e %e ) for 9 G [a,/3]. Of course, we have ip~ 1 (S) = ^ 1 (S-)LlCLIifj^ 1 (S + ) 
as a disjoint union. 

Let {T n (x)} be a sequence of analytic functions on ip~ 1 (S) where we assume that the analytic 
functions satisfy the basic asymptotic relation: 

(17) T n (x) = 1 + a n (x)V>(x) c " + e n (x), 

where {c„} is a increasing unbounded sequence of positive numbers, 5 > is a constant so that 
|an(^)| > 8, and |a n (x)| = exp[o(c n )], uniformly on ip~ 1 (S). The term e n (x) satisfy the following 
estimates uniformly: 

o(ip(x) Cn ), x G S+, 

o{l), xeS-. 

In the sequel, we may assume either form for e n (x) if x lies on the common boundary C of the two 
regions S± that is, \4>{x)\ = 1- 

Let Z n be the set of all zeros of T n that lie in if;~ 1 (S), which we assume is finite for all n. For 

[71,72] C [a, (3), let 

(18) ^(71,72) = #{x G Z n : argx G [71,72]}- 
Choose e > so 3e < eo- By the Argument Principle, we find that 

where T is the boundary of the sector {pe td : p G [1 — e, 1 + e], 9 G [71, 72]}- The closed contour T 
naturally has four parts of the form T\± e and T 7j , j = 1,2 where 

r 1±e = {(l±e)e w :9€ [71,72]}, 
T 7 = {pe*:p€ [l-e,l + e]}. 

A. 2. Contributions over Arcs. 



e n (x) 



„ <2 

lim / 

n ^°° ./rv 



Lemma 9. Inn , rff = . 

Proof. Let £ G ri_ e . Then G V ,_1 (5'-) and e n (x(£)) — ► uniformly on Fi_ e . Furthermore, we 
find 

\a n {x{C)W{x{Q))\ < exp[o(c n )]|C| c " = e~ c " in(i/(i- e )) +o( c„)) = o(1) 
From the basic asymptotic expression 

T n (x(C)) = 1 + a n (x(C))C c " + e n (x(0), 

we find that ln[T n (x(£))] converges uniformly to on ri_ e . Hence their derivatives must also 
converge uniformly to and so the desired integrals converge to 0. □ 
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Lemma 10. lim 9 

n—+ oo 



1 



iT n (x(0) 



72 - 7l ■ 



... T n (x(Q) 

Proof. For ( G Ti+e, |£| = 1 + e and x(Q G V ,_1 (5'+)- By the basic expansion 

T n (x(C)) = 1 + a n (x(C))C c ' 1 + e n (x(C)), 

we have that 



1 + 



l + e n (x(0) 



1 + 



tan 



a n {x{C))C Cn a n (x(C))C c " On(a?(C)) 

We recall our assumptions that |a n (x(C))| > S and e n (x(C))/C Cn = o(l). Hence, we find that 



l + e n (x(0) 



uniformly on T\+ € . In particular, 



r 



+ 



an(x(()) 



Tn(x(Q) 



< 



(1 + e) c " + o(l) 



o(l) 



Tn(x(Q) 



a„(x(C))C e " conver g es to 1 uniformly on T 1+e so M an g ( ' c))C c„ , 
converges uniformly to there as do their derivatives. In other words, we know that 



prMQ) jan(x(Q) 
T„(x(C)) a n (x(C)) 
Note that |a n (x)| = exp[o(c ra )] implies that \a' n (x) / a n (x)\ 

|T n (x(C)) |a n (x(C)) 



^-0 

: o(c n ). Hence, we find 



WO) 



d_ 

dx 



an(x(C)) + f 
a^x)dx + c 
a n (x) c?C C 

T + o(Cn). 



We now conclude that 



i f WO) 
c n r n (x(C)) 



i 



uniformly on Fi +e . 



The lemma now follows easily. 



□ 



A. 3. Backlund's Method. Our estimates for the integrals over T 7 , for 7 G [71,72], are inspired 
by the 1918 method of R. Backlund's proof of the Riemann-von Mangoldt asymptotic formula for 
the number of zeros of the Riemann zeta function. We follow the exposition of Chandrasekharan 
[3] (pages 35-38). 

Recall that T 7 is parametrized as pe iy for p G [1 — e, 1 + e]. 

Let t be the number of zeros of 5ft[T n (»(£))] for ( G T 7 exclusive of endpoints. Then the contour 
integral along T 7 can be written as a sum of integrals over line segments C a b paramerized as pe 11 , 
p G [a, b] where a and b are two consecutive zeros of U[T n (x(0)}- In particular, M[T n (x(pe ie ))] has 
constant sign for p G [a, &]. Then 

|r n (x(C)) 




WO) 



where C n is the image of the line segment pe* 7 , for p G [a, b], under the map £ 1— » T„(x(0)- By 
assumption, C n can only intersect the imaginary axis iR only at its two endpoints; in particular, 
C n must lie either in the left and right half-plane. By Cauchy's theorem, we can deform C n into a 
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semicircle K a b that lies in the same half-plane and has the same endpoints on iM. so that the value 
of the contour integral is unchanged. This allows us to make the estimate 



K ah £ 



< 7T. 



We summarize this discussion as: 
d 

Lemma 11. 



Tn(x(C)) 



p G [1 - e, 1 + e]. 



< (£ + l)vr where t is the number of zeros of ^[T n {x{pe' bS ))\ for 



We now use Jensen's formula to make a useful estimate for £. 



For an analytic function h(z) on some domain E, define h(z) on E c as h(z) = h(z) where 
E c = {z : z G E} which will be analytic on E c . 

For £ G D(l — e; 2e), which is symmetric about the real axis, let 



Tn(0 



T n (x(^))+f n (x(Ze-^)) 



Then f n (Q is analytic on D(l — e; 2e) and 

r„(0 = K[T n (x(£e 17 ))], ee[l-e,l + e]cR. 

For convenience, we recall Jensen's inequality. Let h{z) be an analytic function on the closed disk 
D(a; R), and let < r < R. Suppose h(a) / 0. Let m be the number of zeros of h(z) in the closed 
disk D(a; r) counted according to their multiplicity. Then 



< 



max{ \h(z)\ : \z — a\ = R} 
\Ha)\ ' 



Now each zero of 3?[T n (x(C))] for ( G T 7 corresponds to a zero of f n (Q for £ G [1 — e, 1 + e]. Let £ 
be the number of zeros of T n (£) for £ G D(l — e; 2e). Then we have at once the inequality 

£ < £. 

We will apply Jensen's inequality to the disk D(l — e; 3e) and r = 2e to obtain 



< 



max{|T n (0|:|e- (1-6)| =36} 



|r„(i-e)| 

Lemma 12. (a) f n (l - e) = 1 + o(l). 

ft) max{|f n (0 : |£ - (1 " e)| = 3e} = O (e°^(l + 2e) c ") . 

Proof, (a) Since (1 — e)e 17 G S-, we have the estimate 

|T n (x((l - e)e^))\ < 1 + e°^\l - e) c " + o(l) = 1 + o(l). 

In particular, |T n (l — e)| = 1 + o(l). 

(b) To estimate the maximum of |T n (£)| for |£ — (1 — e)| = 3e, we observe that 

\f n (0\ < \ [\T n (x(^r))\ + \f n (x(te-*))\ . 

Let C he in the closed disk D(l — e; 3e) so \(\ < 1 + 2e. By the basic asymptotic relation, we find 
that whether ( G 5 + or 5_: 

\T n (x(Q)\ < |l + a„(x(C))C c " + e„(x(C))| 

< 1 + e o( - c "\l + 2e) c " + o((l + 2e) c ") = O (e o{Cn) (1 + 2e) c ") , 
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where the big-oh constant holds uniformly on D{1 — e; 3e). In particular, this estimate holds for £ 
that lie on the circle |£ — (1 — e)| = 3e. A similar estimate holds for |T n (x(£e* 7 ))|. We sum up this 
discussion as 

max{|T n (£) : ^ - (1 - e)| = 3e} = O (e°^\l + 2e) c ") . 

□ 

The last lemma together with Jensen's inequality allows us to make an estimate for £: 

( 3 V < ™*{\Tn(0 ■■ \t - (I - e)\ = 3e} = , {Cn) 
V2y " f n (l-e) V 1 ' 

Recalling that £ < £, we have the bound 

ln[0 (e°( c ")(l + 2e) c ™)] 
l ~ M3/2) ' 

Hence, we have the estimate for the integral 



1 



dC 



r 7 T n (x(()) 

Since these bounds hold for all e > sufficiently small, we have shown the 
Lemma 13. For all e > sufficiently small, 

jT n (x(Q 

where T 7 is the line segment pe %1 , p G [1 — e, 1 + e]. 



1 

On 



Finally, this last inequality allows us to make the estimates 
72 " 71 1 ln(l + 2c) + O ( < iV " (7l ' 72) < + -4-- ln(l + 2e) + O ^ 1 



2vr 21n(3/2) v 7 V c n/ ~ 2vrc n ~ 2?r 21n(3/2) v ' \c 

A. 4. Main Density Theorem. By combining the above lemmas and noting that these bounds 
hold for all e > sufficiently small, we obtain our main density result. 

Theorem 6. Let a < 71 < 72 < j3, and let ^(71,72) denote the number of zeros of T n (x) whose 
arguments lie in [71,72], given in equation fi!8\ ). Then 

v ^(71,72) 72-71 
lim = ; 

that is, the image of the zero density under tp is Lebesgue measure on an arc of the unit circle. 

We need to recall the notions of lim sup and liminf of a sequence {X n } of compact sets in the 
complex plane. Now x* S lim sup X n if for every neighborhood U of x, there exists a sequence 
x n , £ X nk n U that converges to x* while x* S lim inf X n if for every neighborhood U of x, there 
exists an index n* and a sequence x n G X n n U, for n > n* that converges to x* . 

It is known that if the liminf X n and lim sup X n agree and are uniformly bounded, then the 
sequence {X n } converges in the Hausdorff metric. 

When the density result holds, then the liminf Z(T n ) must agree with limsupZ(T n ). Hence, we 
have the following: 

Corollary 6. As compact subsets of tp~ 1 (S), Z(T n ) converges to the unimodular curve C in the 
Hausdorff metric. 
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Although we can determine the zero attractor and the zero density completely in the above 
framework, it is conceptually useful to have the result of Sokal that gives a description of the 
support of the zero density measure. 

[Sokal] [6]: Let D be a domain in C, and let zq £ D. Let {g n } be analytic functions on D, 
and let {a n } be positive real constants such that {\g n \ an } are uniformly bounded on the compact 
subsets of D. Suppose that there does not exist a neighborhood V of zq and a function v onV that 
is either harmonic or else identically -co such that lim inf a n In \g n (z) \ < v(z) < lim sup a n In \g n (z) \ 

n >oo n — 

for all z £ V. Then z$ £ lim inf Z{g n ). 

Remark: We can state the asymptotic form for T n {x) in a more symmetric form as: 

N 

T n (x) = V>o0) + '^2a rijk (x)ip k (x) Cn + e n {x) 

k=l 

where N is fixed and the error term has the form 

&n{x) = o(max{^(s) Cn ,0 < k < n}) 

This version explains the asymmetry in our first result where we have ipo{x) = 1 and the zeros 
accumulate along the curve l^o^)! = IVK^OI- 

A. 5. Special Case. Theorem [6] shows that the images of zeros under the conformal map ip are 
uniformly distributed along the corresponding circular arc. It can be applied to many cases that 
arise in a broad spectrum. A special case is worthy of attention; namely, the analytic arc C is a 
straight line segment and i^( x ) ( see Lemma [5J has the form e ax+b , where a and b are constants. 

Corollary 7. If the analytic arc C is a straight line segment and ip(x) is of the form e ax+b , where 
a and b are constants, then the zero density along the line segment C is a multiple of Lebesgue 
measure. 

Proof. Let ci and C2 be the endpoints of the line segment so it is parametrized as x{t) = ci(l— t)+C2t, 
< t < 1. Then tp(x(i)) becomes tp(x(t)) = e <*[ci(i-i)+c 2 t]+6 Since \^( x (t))\ = 1, ip(x(t)) can be 
written as 

V>(x(t)) =e 2ni9 ^, 

where 9{t) is a linear function of t. By Theorem [6] the density of images of zeros under ifj(x) along 
the corresponding circular arc is Lebesgue measure. Hence, its pull-back under ^ is also Lebesgue 
measure since 6{t) is linear. □ 
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